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ABSTRACT 

Context. The measurement of line broadening in cool stars is in general a difficult task. In order to detect slow rotation or weak 
magnetic fields, an accuracy of lkms~' is needed. In this regime the broadening from convective motion become important. We 
present an investigation of the velocity fields in early to late M-type star hydrodynamic models, and we simulate their influence on 
FeH molecular line shapes. The M star model parameters range between log g of 3.0 - 5.0 and effective temperatures of 2500 K and 
4000 K. 

Aims. Our aim is to characterize the T e g- and log g-dependence of the velocity fields and express them in terms of micro- and macro- 
turbulent velocities in the one dimensional sense. We present also a direct comparison between 3D hydrodynamical velocity fields 
and ID turbulent velocities. The velocity fields strongly affect the line shapes of FeH, and it is our goal to give a rough estimate for 
the log g and r cff parameter range in which 3D spectral synthesis is necessary and where ID synthesis suffices. Eventually we want to 
distinguish between the velocity-broadening from convective motion and the rotational- or Zeeman-broadening in M-type stars which 
we are planning to measure. For the latter FeH lines are an important indicator. 

Methods. In order to calculate M-star structure models we employ the 3D radiative-hydrodynamics (RHD) code C0 5 B0LD. The spec- 
tral synthesis on these models is performed with the line synthesis code LINF0R3D. We describe the 3D velocity fields in terms of 
a Gaussian standard deviation and project them onto the line of sight to include geometrical and limb-darkening effects. The micro- 
and macro-turbulent velocities are determined with the "Curve of Growth" method and convolution with a Gaussian velocity profile, 
respectively. To characterize the logg and dependence of FeH lines, the equivalent width, line width, and line depth are regarded. 
Results. The velocity fields in M-stars strongly depend on logg and r cff . They become stronger with decreasing logg and increasing 
r c ff- The projected velocities from the 3D models agree within ~ 100 ms -1 with the ID micro- and macro-turbulent velocities. The 
FeH line quantities systematically depend on logg and T c ff. 

Conclusions. The influence of hydrodynamical velocity fields on line shapes of M-type stars can well be reproduced with ID broad- 
ening methods. FeH lines turn out to provide a mean to measure log g and in M-type stars. Since different FeH lines behave all 
in a similar manner, they provide an ideal measure for rotational and magnetic broadening. 

Key words, line: profiles - stars: low-mass, brown dwarfs - Hydrodynamics - Turbulence 



1. Introduction 

Most of our knowledge about stars comes from spectroscopic 
investigation of atomic or molecular lines. In sun-like and hot- 
ter stars, the strength and shape of atomic spectral lines provides 
information on atmospheric structure, velocity fields, rotation, 
magnetic fields, etc. Measuring the effects of velocity fields on 
the shape of spectral lines requires a spectral resolving power 
between R ~ 10, 000 (Av = 30 km s _1 ) for rapid stellar rotation, 
R <; 30,000 (Av = lOkms -1 ) for slower rotation and high tur- 
bulent velocities, and resolution on the order of R ~ 100, 000 for 
the analysis of Zeeman splitting and line shape variations due to 
slow convective motion. 

In slowly rotating sun-like stars, usually a large number of 
relatively isolated spectral lines are available for the investiga- 
tion of Doppler broadened spectral lines. These lines are embed- 
ded in a clearly visible continuum allowing a detailed analysis 
of individual lines at high precision. At cooler temperature, first 
the number of atomic lines is increasing so that more and more 
lines become blended rendering the investigation of individual 
lines more difficult. At temperatures around 4000 K, molecular 



lines, predominantly VO and TiO, start to become important. At 
optical wavelengths, molecular bands in general consist of many 
lines that are blended so that the absorption mainly appears as 
an absorption band; individual molecular lines are difficult to 
identify. At temperatures in the M type stars regime (4000 K 
and less), atomic lines start to vanish because atoms are mainly 
neutral and higher ionization levels are weakly populated. Only 
alkali lines appear that are strongly affected by pressure broad- 
ening. Thus, the detailed spectroscopic investigation of velocity 
fields in M dwarfs is very difficult at optical wavelengths. 

M-type stars emit the bulk of their flux at infrared wave- 
lengths redward of 1 fim. This implies that observation of high 
SNR spectra in principle is easier in the infrared. Furthermore, M 
type stars exhibit a number of molecular absorption bands in the 
infrared, for example FeH. In these bands, the individual lines 
are relatively well separated and provide a good tracer of stellar 
velocity fields. The lines are intrinsically much narrower than 
atomic lines in sun-like stars because Doppler broadening due 
to the temperature related motion of the atoms and molecules is 
much reduced. Thus, the lines can be used for the whole arse- 
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nal of line profile analysis that has been applied successfully to 
sun-like stars over the last decades. 

Examples of analyses using FeH lines are the investigation 
of the rotation activity connection in field M-dwarfs, which re- 
quires the measuremen t of rotationa l line broadening with an 
accuracy of lkms -1 dReinersI 120071) . Another example is the 
measurement of magnetic fields comparing Zeeman broadening 
in m agnetically sensitive an d insensitive absorption lines (see 
e.g. iReiners & Basrill2006l) . A precise analysis of FeH lines, 
however is only possible if the underlying velocity fields of the 
M dwarfs atmospheres are thoroughly understood. In this paper, 
we model the surface velocity fields of M type stars and their 
influence on the narrow spectral lines of FeH . 

W e calculate 3D-C0 5 BOLD structure models dLudwig et alj 
2002) which serve a s an input for the fin e formation program 
LINF0R3D (based on lBaschek et al.lll966l) . Turbulence's are in- 
cluded in a natural way using hydrodynamics, so that we are able 
to investigate the modeled spectral lines for effects from micro- 
and macro-turbulent velocities in the classical sense and their 
influence on the line shapes. The comparison with ID-models 
gives a rough estimate of the necessity of using 3D-models in 
the spectral domain of cool stars. In the first part of this paper 
we investigate the velocity fields in the models and their depen- 
dence on logg and T e g. In the second part, we investigate the 
influence of velocity fields, log g, and T e ff on the FeH molecular 
lines. 



2. 3D model atmospheres 

The three-dimensional time-dependent model atmospheres 
(hereafter "3D models") are calculated with the radiation- 
hydrodynamics code C0 5 B0LD (abbreviation for "Conservative 
COde for the Computation of COmpressible COnvection in a 
BOx of L Dimensions with L=2,3"). It is designed to model so- 
lar and stellar surface convection. For solar-like stars like the 
M-type objects considered here, C0 5 B0LD employs a local set- 
up in which the governing equations are solved in a small (rela- 
tive to the stellar radius) Cartesian domain located at the stellar 
surface ("box in a star set-up"). The optically thin stellar pho- 
tosphere and the upper-most part of the underlying convective 
envelope are embedded in the computational domain. C0 5 B0LD 
solves the coupled non-linear equations of compressible hydro- 
dynamics in an external gravitational field in three s patial di- 
mensions (Frevtag et al. 2002; Wede mever et al.l 2004) together 
with non-local frequency-dependent radiative transfer. In these 
3D models convection is treated without any assumptions like 
in ID mixing-length theory. The velocity fields and its related 
transport properties are a direct result of the solution of the hy- 
drodynamic equations. Due to this, C0 5 B0LD is a well-suited tool 
to investigate the influence of velocity fields on spectral line 
shapes. A C0 5 B0LD model consists of a sequence of 3D flow 
fields ("snapshots") representing the temporal evolution and spa- 
tial structure of the flow. To perform spectral synthesis calcula- 
tions based on the 3D C0 5 BOLD-models we use the 3D line for- 
mation code Linfor3D. It takes into account the full 3D thermal 
structure and velocity field in the calculation of the line profiles. 
It assumes strict Local Thermodynamic Equilibrium (LTE). In 
this paper we will call the spectral lines computed from three- 
dimensional atmosphere models "3D-lines". 

In order to analyze the influence of velocity fields in M- 
stars on FeH lines, we construct a set of C0 5 BOLD-models with 
r eff = 2500K - 4000K and logg = 3.0 - 5.0 [cgs]. TableQ]give 
the model parameters. In the r e jf-sequence, we simulated main 
sequence stars and varied the surface gravity slightly with in- 



creasing effective temperature. For the log g-sequence, we com- 
puted models with different log g values aiming at the same ef- 
fective temperature of 3300 K but the models settle to slightly 
higher or lower T e g values. We decided not to adjust these re- 
sulting effective temperatures, because slight differences in T e s 
do not change the line profiles significantly. We accepted the T e s 
deviations to avoid the large computational effort which would 
be necessary to adjust the models to a common effective temper- 
ature. However, we apply corrections to the line shape related 
quantities such as equivalent width (see Section FOI l. 

The opacities used in the CO 5 BOLD model calculations 
originate from the PHOE NIX stellar atmosphere package 
( Hauschild t'& Baron|[l999l) assuming a so lar chemical compo- 
sition according to lAsplund et al. ( 20051). The opacity tables 
were c omputed after Fergus on et alj d2005l) and iFrevtag et al.l 
d2009l) . These opacities are particularly well-suited for our in- 
vestigation since they are adapted to very cool stellar atmo- 
spheres. The raw data consist of opacities sampled at 62,890 
wavelength points for a grid of temperatures and gas pres- 
sures. For representing the wavelength dependence of the radi- 
ation field in the C0 5 B0LD models the opacities are re-sampled 
into six wavelen g th groups using the opacity binnin g method 
dNordlundl 119821: iLudwigj 1 19921: iLudwig et alJ [l994). In this 
approach, the frequencies which reach monochromatic optical 
depth unity within a certain depth range of the model atmo- 
sphere, will be grouped into one frequency bin on the basis 
of their opacities. For each investigated atmospheric parameter 
combination the sorting of the wavelengths into groups is per- 
formed according the run of monochromatic optical depth in 
a corresponding PHOENIX ID model atmosphere. The thresh- 
olds for the sorting are chosen in logarithmic Rosseland opti- 
cal depth as {+oo, 0.0,-1.0, -2.0, -3.0, -4.5, -co}. In each group 
a switching is done from a Rosseland average in the optically 
thick regime to a Planck average in the optically thin regime, ex- 
cept for the group representing the largest opacities, where the 
Rosseland average is used throughout. In this last bin, which de- 
scribes the optically thick regions, only the Rosseland average is 
used because the radiati ve transfer is loca l and can be described 
as a diffusive process dVogler et alj|2004l) . 

The horizontal size of the models provide sufficient space to 
allow the development of a small number (10-20) of convective 
cells. Their number has to be large enough to avoid box-size de- 
pendent effects, but also small enough that there is a sufficient 
number of grid points available to resolve each individual cell. 
The size of the convective cells scales roughly inversely propor- 
tional to the surface gravity. Accordingly, the horizontal size of 
the computational box is set to larger sizes towards lower log g 
values. The horizontal size of the model with T e g = 3275 K is 
just large enough to fulfill the criteria of the minimal number 
of 10 convective cells (see Fig. [5]), and we saw in test simu- 
lations that the results will not change with a larger model (in 
horizontal size). Therefore we will use this well evolved model 
as well. The vertical dimension is set to embed the optically thin 
photosphere, and a number of pressure scale heights of the sub- 
photospheric layers below. We deliberately keep the depth of our 
models rather small to avoid problems due to numerical insta- 
bilities analogous to the ones encountered and discussed in our 
previous works on the hydrodyna mics of M-type stellar atmo- 
spheres (|Odwigltin|200l|200l). 

For the comparison with ID models, we spatially average 
the 3D-model over surfaces of equal Rosseland optical depth at 
selected instants in time. We call the obtained se quence of ID 
structu res (3D)-model. We follow the procedure of Steffe rTet al.l 
(fl995l) and average the fourth moment of the temperature and 
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Table 1. Overview of different model quantities for models at constant T^g and different log g (upper part) and at constant log g and 
different T e g (lower part). 



Model code 


Size(x,y,z) [km] 


Grid points (nx,ny,nz) 


H„ [km]" 


z-size [H p ] 


T eS [K] 


l0g£ [Cgs] 


d3t33g30mm00wl 


85000 x 85000 x 58350 


180 x 180 x 150 


2821 


20.7 


3240 


3.0 


d3t33g35mm00wl 


28000 x 28000 x 11500 


180 x 180 x 150 


826 


13.9 


3270 


3.5 


d3t33g40mm00wl 


7750 x 7750 x 1850 


150 x 150 x 150 


250 


7.4 


3315 


4.0 


d3t33g50mm00wl 


300 x 300 x 260 


180 x 180 x 150 


18 


14.5 


3275 


5.0 


d3t40g45mm00n01 


4700 x 4700 x 1150 


140 x 140 x 141 


109 


10.6 


4000 


4.5 


d3t38g49mm00wl 


1900 x 1900 x 420 


140 x 140 x 150 


36 


11.7 


3820 


4.9 


d3t35g50mmOOwl 


1070 x 1070 x 290 


180 x 180 x 150 


20 


14.5 


3380 


5.0 


d3t28g50mm00wl 


370 x 370 x 270 


250 x 250 x 140 


13 


20.8 


2800 


5.0 


d3t25g50mm00wl 


240 x 240 x 170 


250 x 250 x 120 


12 


14.2 


2575 


5.0 


" at r = 1 



first moment of the gas pressure to preserve the radiative proper- 
ties of the 3D-model as far as possible. The 3D velocity informa- 
tion is ignored in the (3D)-model and replaced by a micro- and 
macro-turbulent velocity. By construction, the (3D)-model has 
the same thermal profile as the 3D-model, but evidently without 
the horizontal inhomogeneities related to the convective granu- 
lation pattern. We will call the spectral lines synthesized from 
<3D)-models "<3D>-lines". 

2.1. Atmosphere structures 

The temperature stratification shown in FigQ] (top left) of the 
models with changing T e s appears very similar for all models 
in the region below log t ~ I. They reach their T e s around 
t ~ 2/3 and continue to decrease to higher atmospheric layers. 
Above log r ~ 1 the temperature of the two hottest models in- 
crease more strongly than in the cooler cases. Since the models 
are almost adiabatic in the deeper layers (see below) the tem- 
perature gradient follows the adiabatic gradient, which is given 
by the equation of state and steeper in hotter models due to the 
inefficient H2 molecule formation. This increase of temperature 
is also very similar to the lower log g models on the right side 
of FigQ] These models also show a temperature gradient which 
becomes steeper to deeper atmospheric layers and with decreas- 
ing surface gravity again due to a steeper adiabatic gradient. All 
models reach an almost equal effective temperature at t ~ 2/3 
and their temperature stratification does not differ strongly to 
smaller optical depth. 

The entropy stratification for models with varying (mid 
left in FigQ]) show a similar behavior for all models. It is adia- 
batic (dS jdx — 0) in layers below logr ~ 1 and has a supera- 
diabatic region (dS/dr > 0) between logr ~ 1 and ~ - 1 
that moves slightly to smaller optical depth for hotter models. 
In these regions, with dS /dr > 0, the models are convectively 
unstable and become convectively stable in the outer parts of 
the atmosphere where dS/dr < 0. In the models with chang- 
ing logg (mid right in FigQ}, we can see that the entropy be- 
haves almost like in the T e ff case. To lower surface gravities, 
the superadiabatic region is more significant. In higher layers, 
the models become convectively stable except the model with 
logg = 4.0 [cgs] which shows a second decrease of entropy to 
the outer part layers. To understand this behavior, we have to 
investigate the adiabatic gradient of this region which is very 
small and changes very little along the upper atmosphere. This 
is due to the equat ion of state used i n the models and can be 
seen in Fig. 16 of iLudwig et al.l ([2006) (model H4 in this figure 
equates to our logg — 4.0 model). This figure shows that the 
upper atmosphere lies in a plane of small and constant adiabatic 



gradient. Due to this, the model becomes convectively unstable 
again in the upper layers. This is probably the reason for the, in 
comparison to other models, higher velocities in this model. 

In the left bottom panel of FigQ] the horizontal and vertical 
rms-velocities are plotted for models with different T e ff . Both ve- 
locity components increase with increasing T e ff . The maxima of 
the vertical velocity moves to slightly deeper layers with higher 
temperatures and the maxima of the horizontal velocity stays al- 
most at the same optical depth. We can see a qualitatively similar 
dependence in the log g model sequence in the right bottom plot 
in FigQ] Only the model with logg = 4.0 [cgs] shows peculiar 
behavior in the upper atmospheric layers, which is probably re- 
lated to the entropy stratification in this model. We will describe 
the velocity fields in the models in more detail and with a slightly 
different method in the next section. 

3. Velocity fields in the C0 5 BOLD-models 

Before we investigate the effect of velocity fields on spectral 
lines, we analyze the velocity fields in the models themselves 
and we will do this with respect to spectral lines. Spectral lines 
are broadened by velocity fields where the wavelength of absorp- 
tion or emission of a particle is shifted due to its motion in the 
gas. Here we are mostly concerned with the macroscopic, hydro- 
dynamic motions but have in mind that the thermal motions are 
also constituting an important contribution. If we envision each 
voxel in the RHD model cube to form its own spectral line, the 
whole line consists of a (weighted) sum of single lines. The ve- 
locity distribution might be represented by a histogram of the ve- 
locities of the voxels which gives us the velocity dispersion. We 
try to describe the velocity fields in that sense instead of using 
the rms-velocities shown in FigQ] In the C0 5 BOLD-models, a ve- 
locity vector is assigned to each voxel and consists of the veloci- 
ties in x-, y-, and z-direction. We will investigate the vertical and 
horizontal component of the velocity dispersion in the models 

and the total velocity dispersion cr tot = Jcr\ + cr^ + cr\. In order 

to describe the height dependent velocity dispersion we applied a 
binning method, i.e. we plot all velocity components of a certain 
horizontal plane of equal optical depth r in the C0 5 B0LD cube in 
a histogram with a bin size of 25 m s _1 . We fit the histogram ve- 
locity distribution with a Gaussian normal distribution function 
G = exp(-(^) 2 ) and take the standard deviation cr as a measure 
for the velocity dispersion cr in the models (see Fig. |2j (The rela- 
tion between the Gaussian standard deviation cr and the standard 
deviation <x rms of the mean velocity is cr — V2 ■ cr rms ). This is 
done for the cr x , <x v , and cr, component of the velocity vector for 
each horizontal plane from Tmi n to T max , which are the highest 
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Fig. 3. Upper panel: Radial (o%) and horizontal (cr H ) components of the velocity dispersion are plotted against the optical depth 
on a logarithmic scale. Bottom panel: Each bottom panel shows the contribution functions (CF) of an FeH-line at a wavelength of 
9956.72 A. Equivalent width W,\ (solid) and the depression at the line center (dashed) of the line are plotted as a function of optical 
depth on a logarithmic scale. The models (from left to right) are located at of 2800 K, 3380 K, 3820 K and 4000 K and a logg 
value of 5.0, except the one with T eS = 3820K(logg = 4.9), and the one with T eff = 4000K(logg = 4.5) [cgs]. 
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Fig. 4. Upper panel: Radial (o%) and horizontal (cr H ) components of the velocity dispersion are plotted against the optical depth 
on a logarithmic scale. Bottom panel: Each bottom panel shows the contribution functions (CF) of an FeH-line at a wavelength of 
9956.72 A. Equivalent width W,\ (solid) and the depression at the line center (dashed) of the line are plotted as a function of optical 
depth on a logarithmic scale. The models are located at log g values from left to right of 5.0 (T e ff = 3275 K), 4.0 (T^ = 3315 K), 
3.5 (T eS = 3270 K) and 3.0 (T eS = 3240 K) [cgs]. 



and the deepest points, in each model atmosphere. In this way 
we obtained the height dependent velocity dispersion cr x y Jj). 
We average over five model snapshots for a better statistical sig- 
nificance. In Figs.[3]and|4] the velocity dispersions for the hori- 
zontal components cth = Jcr* + &y and vertical component cr z 

are plotted against optical depth (black solid and dotted lines). 
In the latter figures, we identify the maxima around log t = — 1 
of the velocity dispersions as the region where the up-flowing 
motion spreads out in horizontal direction and starts to fall back 
to deeper layers. We will call this point, the "convective turn- 
over point". In a horizontal 2D cross-section of the vertical ve- 
locities, somewhat below this area, the up-flowing granulation 
patterns are clearly visible (see Fig|5]for models with different 



log g or Ludw iget al.l d2006l) ). In a vertical 2D cross-section one 
can see the behavior of the up-streaming material (Fig.|6]shows 
an example for a model with 7^ = 2800 K). The material starts 
moving upwards almost coherently and before reaching higher 
layers (around z - 200 km in Fig. |6]l of the atmosphere the ver- 
tical velocity dispersion cr z becomes maximal. After that point, 
the material spreads out in horizontal directions and starts falling 
down again. At this point, the dispersion of horizontal velocities 
<th becomes maximal (convective turn-over point). 

In Fig|4]we can see at lower surface gravities that the max- 
ima of the horizontal velocity dispersion are not centered around 
a specified optical depth anymore; it spreads out in vertical direc- 
tion and span the widest range at log g — 3 .0 [cgs] . The pressure 
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Fig. 1. From top to bottom, the temperature, entropy, vertical, 
and horizontal velocity are plotted as a function of optical depth. 
The column on the left side shows models with different T e g 
and constant logg, on the right side the models are at a con- 
stant T e ff with varying logg. The rise of horizontal velocity in 
the log g = 4.0 [cgs] model in the deeper atmospheric layers is 
due to interpolation from z- to r-scale. 
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Fig. 5. Horizontal cross-section around t ~ 1 of vertical veloc- 
ity components. The models are located at \ogg values of 3.0, 
3.5, 4.0, and 5.0 [cgs] (from upper left corner to lower right, re- 
spectively). 





velocity in z direction [m/s] 



Fig. 2. Histogram of the velocity distribution in vertical direc- 
tion. The normalized number of points is plotted against the 
vertical velocity in m/s (solid line). The Gaussian (dashed line) 
fits the velocity distribution and determined an FWHM value 
(dashed-dotted line), which is related to <x with FWHM = 
2 Vln2 ■ <x. The underlying model is located at T e g = 2800 K 
and log g = 5 [cgs]. 




u 



Fig. 6. 2D cross-section in x-z direction of the x-z velocity in a 
model with log g = 5.0 [cgs] and T e $ - 2800 K. Overplotted 
are x-z velocity vectors. The direction of the material flow is 
indicated by the velocity vectors. 



stratification changes, and the convective turn-over point moves 
to lower gas-pressure (not shown here) but stays at almost con- 
stant optical depth between logr = and logr = - 2. 
With varying temperature, the position of the convective turn- 
over point stays at almost constant optical depth. 

3. 1 . Reduction of the 3D velocity fields 

Commonly, micro- and macro-turbulence derived from spec- 
troscopy are interpreted as being associated with actual veloc- 
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Fig. 7. Plotted are the weighted projected and unprojected veloc- 
ity dispersions of the horizontal and vertical component for mod- 
els with different T e g. The models are located at T e g values of 
2575 K, 2800 K, 3275 K, 3380 K, 3820 K and 4000 K and a log g 
value of 5.0, except the one with T e ff = 3820 K (log g = 4.9), 
and the one with T eS = 4000 K (log g = 4.5) [cgs]. 



troscopy (see appendix). When interpreting comparisons shown 
below the very approximate nature of our model should be kept 
in mind. In this model, we include the geometrical projection 
of the components of cr XM to the line of sight of the observer. 
We also have to consider the effect of limb-darkening of the 
stellar disk. For each velocity dispersion component we calcu- 
late a projection factor which includes both, geometrical pro- 
jection and limb-darkening effects (more described in the ap- 
pendix). We take a limb-darkening coefficient of 0.4 which fol- 
lows from the continuum from angle dependent line synthesis 
performed in LINF0R3D. These simulations suggest that a linear 
limb-darkening law with a limb-darkening coefficient of 0.4 is 
suited to describe the brightness variation. The projected veloc- 
ity dispersions are also plotted in Figs.[3]and[4](red solid and dot- 
ted lines). The reducing effect of this projection factor is stronger 
in the horizontal components than in the vertical because the pro- 
jected area at the limb of the stellar disk, where cr H reaches its 
maximum value, is much smaller than in the center where cr z 
has its maximum value. The influence of limb-darkening is not 
strong and the dependence of the projection factor from the limb- 
darkening coefficient is only small (described in more detail in 
the appendix). 



> 1000- 
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Fig. 8. Plotted are the weighted projected and unprojected veloc- 
ity dispersions of the horizontal and vertical component for mod- 
els with different logg. The models are located at T e s around 
3300/T and different log g values of 3.0 (T eS = 3240 K), 3.5 
(r eff = 3270 K), 4.0 (T eS = 3315 K) and 5.0 (r eff = 3275 K) 
[cgs]. 



ity fields present in the stellar atmosphere. In our simulations no 
oscillations are induced externally but small oscillations are gen- 
erated in the simulations itself. The velocity amplitudes of these 
oscillations reach maximal 10 % of the convective velocities and 
have no significant influence on the macro-turbulent velocity. We 
would also not expect global oscillations for these objects, ex- 
cept for young stars with solar masses smaller than 0.1 Mq 
induced by D-burning dPalla & Baraffd 120051) . In the follow- 
ing, we try to make the connection between micro- and macro- 
turbulence and actual hydrodynamical velocity fields manifest 
by considering the velocity dispersion determined directly from 
the hydrodynamical model data, and compare it with the micro- 
and macro-turbulence derived from synthesized spectral lines 
(see Sec. 3. 2). This connection is algebraically not immediate, 
and we only apply a simple model for translating the hydro- 
dynamical velocities into turbulent velocities relevant to spec- 



3.1 .1 . Weighted velocities 

To investigate the influence of broadening from the projected 
and unprojected velocity dispersion on spectral lines, we use 
contribution functions for the equivalent width Wa and the de- 
pression at the line center of an FeH line at 9956.7 A (Magain 
1986). The line gains its Wa and depression in the region be- 
tween logr = 1.0 and logr = - 4.0, i.e. that is the re- 
gion of main continuum absorption caused by FeH molecules. 
The maximum is roughly centered around log r = -1.0 and 
moves to slightly lower optical depth with lower temperatures 
(at the lowest T e g of 2575 K, the maximum is centered around 
logr = - 2.0) or higher surface gravities. The contribution 
function of Wa range over the region of the convection zone and 
reflects its influence on the line shape. Due to the latter fact, 
FeH lines are a good mean to explore the convective regions in 
M-dwarfs. In order to measure the velocities in the region where 
the lines originate, we compute the mean of the (projected) ve- 
locities weighted with the contribution function of Wa- 



C„ e ighted - 6 
2j t =2 



(1) 



The horizontal and vertical components of these weighted ve- 
locity dispersions are plotted in Figs. [7] and [8] We can see an 
increase of cth,- with increasing effective temperature or decreas- 
ing surface gravity. We can see again that the projected horizon- 
tal velocity dispersions are significantly smaller than the unpro- 
jected ones due to the reasons mentioned above. The difference 
in the vertical component is much smaller. The velocity disper- 
sions for the vertical component range from a few hundred m/s 
for cool, high gravity models to ~ 1 km s _1 for hot or low grav- 
ity models. The horizontal component range from a ~ 500 m s _1 
for cool, high gravity models to ~ 2 km s _1 for hot or low grav- 
ity models. We compare the total projected velocity dispersion 

Utot - -tJcIi + °"? m Sec. 3.2 with micro- and macro-turbulent 
velocities in the classical sense. 

The strong increase of the velocity dispersions in the atmo- 
spheres to higher layers (Figs. [3] and |4j, which some models 
show, is related to convective overshoot into formally stable lay- 
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ers. These velocities are originated by waves excit ed by stochas- 
tical fluid motions and by advective motions dLudwig et alj 
120021 and references therein). However, it will not affect the 
spectral lines, because the lines are generated in the region be- 
tween an optical depth of log 7 = 1.0 and log r = - 4.0. The 
lines on the model with T e g = 2575 K are an exception, they are 
formed in the outermost layers of the model and it is not possible 
to compute the full range of formation of these lines, because the 
atmosphere is not extended enough. One has to keep this in mind 
when regarding the line dependent results of this model later in 
this chapter. 

3.2. Radial velocity shifts 

Due to the fact that in a convective motion the up-flowing area is 
larger, because it is hotter and less dense, than the down-flowing 
one, one expect a net shift of the velocity distribution to positive 
velocities. That means, the net amount of up-flowing area with 
hotter temperature, i.e. more flux in comparison to the down- 
flowing area, results in a blue shift in the re st wavelength posi- 
tion of a spectral line (see e.g. lDravrnslll982l) . 

To see how the area of up-flowing material affects the rest 
wavelength position of a spectral line, we computed ten FeH 
spectral lines (more described in Sect. 4) on 3D models to mea- 
sure the displacement of the line positions. In order to deter- 
mine the center of the line, we used the weighted mean C = 
which account for the asymmetric line shape. (To use the 
weighted mean is appropriate here, since we have no noise in 
the computed data.) The line shifts of the flux and the inten- 
sity are given in Tab. [2] A negative value stands for a blue shift, 
and a positive for a red shift. The values for each model are the 
mean of five temporal snapshots. The absolute displacement of 
the flux and intensity in the log g series reflects the dependence 
of the velocity fields on surface gravity, but for the T e « series a 
connection is barely visible (see Figs. [7] and [8). This could be 
due to the small geometrical size of the atmospheres in the T e ff 
series (see Tab.[TJ. Only the the one with T e g = 4000 K shows a 
significant line shift and in this model the atmosphere is 1 150km 
high due to the slightly smaller logg value of 4.5 [cgs]. At this 
point we will not go in deeper analysis of this topic. 

Since we only use five snapshots, we are dealing with statis- 
tics of small numbers and hence a big scatter in the results. This 
scatter cr s hift is in general one order lower than the shift of the 
line and the integrated jitter of the line, which is important for 
radial-velocity measurements, goes with (Tj; tter = where N 
is the number of snapshots. Since in a star N is of the order of 
10 6 , the jitter will be of the order of mm s _1 . 

We did not further investigate the effect of granulation pat- 
terns on the line profiles but, as we will see below, the lines 
are almost Gaussian and show no direct evidence for significant 
granulation effects. 

3.3. Micro- and macro-turbulent velocities 

Due to the large amount of CPU-time required to compute 3D 
RHD models and spectral lines on these models, we want study 
the necessity of 3D models in the range of M-stars. Our goal 
is to compare the broadening effects of the 3D velocity fields 
on the shape of spectral lines with the broadening in ter ms of 
the classical micro- and macro-turbulence profiles (see e.g. lGravl 
Il977ll2008h . The latter descripti on is common ly used in ID at- 
mosphere models like ATLAS9 dKuruczlll970T) or PHOENIX and 
related line formation codes. If the difference between ID and 



3D velocity broadening is small, the usage of fast ID atmosphere 
codes to simulate M-stars for comparisons with observations, 
e.g. to determine rotational- or Zeeman-broadening, would be 
an advantage. 

If the size of a turbulent element is small compared to unit 
optical depth, we are in the regime of the micro-turbulence. The 
micro-turbulent velocities might differ strongly from one posi- 
tion to another and have a statistical nature. The broadening ef- 
fect on spectral lines can be describ ed with a G aussian which en- 
ters the line absorption coefficient (Grav 2008). It can be treated 
similar to the thermal Doppler broadening. The effect on the 
shapes of saturated lines is an enhancement of line wings due 
to the fact that at higher velocities the absorption cross section 
increases and as a consequence the equivalent width (Wa) of the 
line is increased. 

If the size of a turbulent element is large compared to unit op- 
tical depth (maybe better of the same size as unit optical depth), 
we are in the regime of the macro-turbulence. It can be treated 
similar to rotational broadening as a global broadening of spec- 
tral lines. The effect is an increase of the line width but the equiv- 
alent width remains constant. 

As we saw before, the velocity fields in the M-stars are not 
very strong in comparison to the sound speed (see Tab|2]i and one 
could expect that their influence on line shapes do not deviate 
strongly from Gaussian broadening. 

We com pared l ine broadening with the radial-tangential 
profile from |GrayJ d 1975b and a simple Gaussian profile and 
found that the latter is a good approximation with an accu- 
racy high enough for possible determination of rotational- or 
Zeeman-broadening. Hence, in this investigation we will assume 
Gaussian broadening profiles. That means, we can assume a 
height independent isotropic velocity distribution for micro- and 
macro-turbulent velocities. This is a very convenient way to sim- 
ulate the velocity fields. One would expect that the anisotropic 
nature and the height dependence of the hydrodynamical veloc- 
ity fields have a significant influence on line shapes, so it is re- 
markable that their influence on spectral lines can be described to 
a high accuracy in this way. (At least in the investigated M-type 
stars.) In Fig. [9] a few examples of (3D)-FeH spectral lines were 
plotted, which were computed with a given micro-turbulent ve- 
locity (determined below) and then convolved with a Gaussian 
broadening profile with a given macro-turbulent velocity. The 
broadened (3D)-FeH lines fit the 3D-FeH lines very well. The 
difference of the ID and 3D centroid (C = range in the 

order of m/s for small velocity fields to 30 - 40 m s _1 for strong 
velocity fields in hot M star models or with low log g. The error 
in flux is lower than 1% (see Fig. [9j» this corresponds to an un- 
certainty in velocity, for example rotational velocity, of less than 
150 ms -1 depending on the position in the line. It is also visible 
in the latter figures that at low effective temperature effects from 
velocity broadening are not visible in comparison with an un- 
broadened (3D)-line in which the v.d.Waals broadening is dom- 
inant. At higher effective temperatures, the difference between 
broadened and unbroadened (3D)-lines is clearly visible. We 
found that in the range of M-type stars, ID spectral synthesis 
of FeH-lines using micro- and macro-turbulent velocities in the 
classical description is absolutely sufficient to include the effects 
of the velocity fields. In the following we will determine the ve- 
locities needed. 



8 S.Wende et al.: 3D simulations of M star atmosphere velocities and their influence on molecular FeH lines 

Table 2. Displacement Avfi ux and Avi ntens ity in m/s of the position of an FeH line from the rest wavelength on models with different 
r e ff (left side) and different logg (right side). 



Teff 


Avpiux [m/s] 


AVjntensity [m/s] 


Vad [m/s] 


l°gg 


Av F i ux [m/s] 


AVintensity [m/s] 


Vad [m/s] 


2800 





2 


3750 


3.0 


-56 


-25 


6100 


3275 


-1 


3 


4000 


3.5 


-47 


-36 


5500 


3380 


-2 


2 


4850 


4.0 


-16 


-7 


5400 


3820 


-1 


-5 


6000 


5.0 


-1 


3 


4500 


4000 


-44 


-28 


7000 












5 
Doppler Velocity [km/s] 




-2 2 

Doppler Velocity [km/s] 

Fig. 9. FeH lines for models with r efF = 2800 K, log g = 5.0 
[cgs] (top) and 7/ eff = 3820 K, log g = 4.9 [cgs] (bottom). The 
upper panels show the 3D-line (dots) and the (3D) convo i ve d-line 
(solid line) which was broadened by a Gaussian profile. For com- 
parison we plotted a (3D)-line which was not broadened by any 
velocities (dashed line). In the lower panels the 3D-(3D) convo i ve( j 
residuals are plotted. One can see the asymmetry which stems 
from the line shifts due to convective motions. Note the different 
y-axis scale. 



3.3.1. Determination of micro- and macro-turbulent velocities 

The investigation of the micro-turbulent v elocities w as done 
with the curve of growth (CoG) method (e.g., Gray 2008). We ar- 
tificially increase the line strength of an absorption-line (increase 



the log gf value) which in turn increase the saturation of the line 
and its influence of the micro-turbulent velocity, which results in 
an enhancement of Wa- In order to determine micro-turbulent- 
velocities, we use Fel- and FeH-lines produced on (3D)-models 
with different micro-turbulent velocities (there are no differences 
in micro-turbulent velocities between both types of lines), i.e. for 
each log g/-value we compute a (3D)-line with micro-turbulent 
velocities between 0kms~' and lkms -1 in0.125kms~' steps. 
In this way we obtained CoGs for 9 different micro-turbulent ve- 
locities. We compare the equivalent widths in the CoGs with the 
ones which was computed on the 3D-models and selected the 
velocity of the CoG which fits the 3D CoG best in the sense of 
X 2 -residuals. The macro-turbulent velocities are determined by 
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t .. run 
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Fig. 10. Macro- (solid) and micro- (dotted) turbulent velocities 
and the sum of both (long dashed line) as a function of T e ff. The 
data points are fitted by a second order polynomial and the sum 
of micro- and macro-turbulent velocities is given by the sum of 
the fits. For comparison, velocities from Bean et al. and the to- 
tal projected weighted 3D velocity dispersions are plotted, too 
(short dashed line). The models are located at T e g values of 
2800 K, 3275 K, 3380 K, 3820 K and 4000 K and a logg value 
of 5.0, except the one with T eff = 3820 K (logg = 4.9), and 
the one with r etf = 4000K(logg = 4.5) [cgs]. 



computing a (3D)-model FeH line (at A = 9956.7 A) including 
the micro-turbulent velocities. This line is then broadened with 
a Gaussian profile and different velocities until it matches the 
3D-model line. 

The dependence of the micro- (£) and macro- (£) turbulent 
velocities on surface gravity and effective temperature is plotted 
in Figs. [TUl and [TTI The macro-turbulent velocities in both cases 
show a quadratic dependence, and we can fit them with a second 
order polynomial. The micro-turbulent velocities could be fitted 
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log(g) [cgs] 
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Fig. 11. Macro- (solid) and micro- (dotted) turbulent velocities 
and the sum of both (long dashed line) as a function of log g. The 
data points are fitted by a second order polynomial and the sum 
of micro- and macro-turbulent velocities is given by the sum of 
the fits. For comparison, the total projected weighted 3D velocity 
dispersions is plotted, too (short dashed line). The models are 
located at 73 e ff around 3300 K and different logg values of 3.0 
(T eS = 3240 K), 3.5 (7/ eff = 3270 K), 4.0 (r eff = 3315 K) and 
5.0(r eff = 3275 K) [cgs]. 



by a linear function or a second order polynomial. We decided 
to use the second order polynomial, too. 

Micro- and macro-turbulence velocities both show a similar 
dependence on surface gravity and effective temperature, which 
implies that there is a direct connection between both. A compar- 
ison of the macro (£)- and micro (£)-turbulent velocities with the 
sum of both ( ^ 2 + £ 2 ) and the total projected weighted veloc- 
ity dispersion cr tot (see Sec. 3) is also shown in Figs.[T0land[TTl 
The total projected weighted velocity dispersion (see Sec. 3) is 
very similar to the macro-turbulent velocities and in most cases 
smaller than the sum of micro- and macro-turbulent velocity. It 
is remarkably, that it is possible, with this simple description of 
the total projected weighted velocity, to describe the broadening 
influence of the hydrodynamical velocity fields on spectral lines 
in comparison with the classical micro- and macro-turbulent de- 
scription. 

In order to obtain a good estimate of the line profile in ID 
spectral line synthesis, knowledge of the micro- and macro- 
turbulent velocities is very important. Otherwise one could un- 
derestimate the equivalent width or the line width and hence ob- 
tain a wrong line depth. 

We compare our micro- an d macro-turbulent ve locities to ob - 
servational results from (see Be an et al. l l2006allbHBear3l2007h . 
Our value of the macro-turbulent velocities are roughly of the 
same order. The higher macro-turbulent velocities from Bean 
et al. possibly contain rotational broadening, but the Bean et 
al. micro-turbulent velocities are roughly a factor two or three 
higher than ours. These velocities are obtained from observed 
spectra by the authors of the afore-mentioned papers using 
spectral fitting procedures. In detail, they used PHOEN IX at- 
mosp here models and the stellar analysis code MOOG dSnedenl 
1 1973b . One has to keep in mind that the empirical determina- 
tion of micro-turbulence is also endangered to suffer from sys- 
tematic errors. For most o f the lines that Bean and collabora- 
tors employ (line data from lBarklem et ai1 d2000)). the v.d.Waals 
damping constant is available. However, in some cases not, 



then usually Unsold's hydrogenic approximation is applied to 
calculate the value, and different authors use significantly dif- 
ferent enhancement factors changing its value. This illustrates 
the level of uncertainty inherent to this approach. For instance, 
ISchweitzer et al.l ([1996) use an enhance factor of 5.3 (for the re- 
sulting y6 valu es) for thei r Fe I lines, while Bean et al. prefer 2.5 
for Til lines (Bean 2007). To investigate the detailed influence 
of v.d.Waals broadening on determination of micro-turbulence 
velocities is very interesting, but goes beyond the scope of this 
paper. We want to point out, that it is perceivable that uncertain- 
ties in the damping constant can introduce significant system- 
atic biases in the resulting value of spectroscopically micro- and 
macro-turbulence and they could easily be overestimated. 

As mentioned above and illustrated Fig. [10] our predic- 
tion of the micro-turbulence grossly underestimates the micro- 
turbulence values measured by Bean et al.. This might hint at 
deficits in the hydrodynamical modeling and we cannot exclude 
the possibility that a process is missing in our 3D models lead- 
ing to a substantially higher micro-turbulence. But due to the ar- 
gumentation above and a comparison with the solar micro- and 
macro-turbulence, we argue that the Bean et al. values for the 
micro-turbulence are too high. However, before being able to 
draw definite conclusions the observational basis has to be en- 
larged. 



4. r eff - and log ^-dependence of FeH molecular lines 

In this section we study the dependence of FeH molecular lines 
on T e g and \ogg in our 3D- and (3D)-models. Again, the inten- 
tion is to identify multi-D effects which might hamper the use 
of the FeH line diagnostics in standard ID analyses. For this 
purpose we compute the (3D)-lines with no micro- and macro- 
turbulence velocity. With this method we can study the FeH lines 
without any velocity effects and can through direct comparison 
between 3D- and (3D)-lines clearly identify velocity induced ef- 
fects. 



4.1. FeH line data 



Win g & Fordl ([1969) were the first to detect a broad molecular 
absorption band arou nd 99 l nm in late M dwarfs. T his band was 
later f ound in S stars dWing||1972l) and in sun spots. Nordh et al] 
dl977l) identified the Wing-Ford band as the 0-0 band of a FeH 
electronic transition. The FeH molecule is well suited for the 
measurements mentioned in the introduction because of its in- 
trinsically narrow and well isolated spectral lines. These lines 
are also an ideal tracer of line broadening in M-stars due to con- 
vection or very slow rotation (Reiners 2007). Since FeH lines 
were not very commonly used for the interpretation of stellar 
spectra in the pa st, only little data is recorded in the literature. 
With the work of Dulick etaLl d2003l) . it is possible to determine 
the log gf value and the transition energies. For the partition 
functio n, a combination of the tabulated function in lDulick et al.l 
(I2003I) and an analyti cally determined one from the Eq.l in 
ISauval & Tatuml dl984l) is used. 

Due to the high atmospheric pressures, v.d.Waals broadening 
is often significant in cool M-type dwarfs. No detailed calcula- 
tions exist for the v.d.Waals broadening of FeH molecular lines. 
Lacking a more accurate trea t ment, we follow the approximate 
approach of ISchweitzer et ail (I1996I) and apply Unsold's hydro- 
genic approximation, albeit - different from Schweitzer and co- 
workers - we do not apply an enhancement factor to the calcu- 
lated broadening constant C(,- 
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The ionization energy of the FeH molecule enters the cal- 
culation of Ce which, unfortunately, is not known. Only the 
dissociation energy is published of 1.59 + 0.08 eV at OK (e.g. 
iDulick et al.l [2003). To derive an estimate of the ionization 
energy, we heuristically compare the ionization and dissocia- 
tion energies of a large number of hydrides taking data from 
Wilkinson] d 19631) . We find an approximately linear relationship 
between the ionization and the dissociation energies of hydrides. 
For FeH we obtain an ionization energy of 6 e V for the known 
dissociation energy from a linear fit. Computation of this value 
from the ionization-potential of Fe and the dissociation energies 
of FeH and FeH + yields 7.3 eV (Bernath 2008, private com- 
munication), which is compatible with our value considering 
our rather crude procedure. The difference in FWHM of syn- 
thesized FeH-lines in the range between 6 and 7.3 eV amounts 
to ss 25 m/s. This uncertainty is quite acceptable in comparison 
to the total broadening of typically several 100 ms 1 in this in- 
vestigation. 



4.2. An ensemble of 3D- and <3D)-FeH lines 

We i nvestigate ten FeH lines between 9950 A and 9990 A chosen 
from Reiners & Basril d2006h (see Tab. [3]). We choose lines from 
different branches (Br), orbital angular momentum u>, and rota- 
tional quantum number J. The wavelengths in Tab.[3]are given in 
vacuum and E/ is the lower transition energy. While not directly 
relevant in the present context, because we do not study the ef- 
fects of magnetic fields, we note that five lines are magnetically 
sensitive and five insensitive. We performed the line synthesis 
for fixed abundances on the C0 5 B0LD atmosphere models listed 
in Tab. Q] The spectral resolution is R « 2 • 10 6 (= 5 ■ 10~ 3 
A) corresponding to a Doppler velocity of v ~ 150 m s _1 at the 
wavelength of the considered lines (~ 9950 A). Figures [12] and 
[13] illustrate the strong influence of surface gravity and effective 
temperature on the line shape for the 3D-models. In Fig. [12] one 
can see that for both, (3D)- and 3D-lines, the line depth, line 
width and equivalent width Wa decrease strongly with increas- 
ing effective temperature. The decrease of Wa is due to stronger 
dissociation of the FeH molecules at higher temperatures, i.e the 
number of FeH molecule absorber decreases. Differences in the 
line shape between 3D- and (3D)-lines with changing temper- 
ature are hardly visible. To higher T e s values, the influence of 
broadening on the 3D lines due to velocity fields is slightly vis- 
ible and not covered by thermal and v.d.Waals broadening any- 
more. To cooler temperatures, the velocity fields decrease and 
the differences between the (3D)- and 3D-line shapes vanishes. 
The v.d.Waals broadening is larger then the thermal broaden- 
ing or that from the small velocity fields in the RHD models. 
In the model with T e g = 2800 K, the lines start to become 
saturated. The FeH lines in the z-band at effective temperatures 
below ~ 2600 K become too saturated and too broad for inves- 
tigations of quantities like magnetic field strength or rotational 
broadening below lOkms -1 . 

The differences between 3D- and (3D)-line shapes of the 
log ^-series in Fig.Q~3]is more obvious than in the r e ff case. The 
differences in line depth and line width become significant with 
decreasing logg. The lines in the 3D-models are significantly 
broadened due to the velocity fields in the RHD models, hence 
the line width is larger and the line depth lower. As we saw in 
Chapter 3, these velocity fields increase with decreasing \ogg 
and could be described in the ID case in terms of macro- and 
micro-turbulent velocities. 



The (3D)-lines become slightly shallower and narrower to- 
wards smaller log g. The equivalent width of the lines decreases 
with decreasing log g due to decreasing pressure and hence de- 
creasing concentration of of FeH molecules. Also the v.d.Waals 
broadening looses its influence to lower pressures and the line 
width decreases. We point out, that we use in all models the same 
solar like chemical compositions, hence the concentration of Fe 
and H stays the same. The creation of FeH also depends on the 
number of H2-molecules, which become larger towards lower 
temperatures and will be important in cool models. 

The slightly different effective temperatures in the models 
with different log g (see Tab. [D affect the line depths as well. 
If the effective temperatures would be the same in the \ogg- 
models, one could expect a monotonic behavior in decreasing 
line strength for decreasing surface gravity in the (3D)-lines. 
However, because the model with logg = 3.0 is cooler, the 
line depth is deeper than that of the one with log g = 3.5. In the 
following analysis we will correct the FWHM, Wa, and the line 
depth of the lines on models with different log g for their slightly 
different effective temperatures. 

The ten FeH lines behave all in the same way as the pre- 
sented ones. We do not see any effect of different excitation po- 
tentials or log gf values in the line shapes that cannot be ex- 
plained by their different height of formation. Thus, we expect 
that we can exclude an extraordinary interaction between these 
quantities and effective temperature or surface gravity. We will 
quantify this preliminary result in the next section. 



Tabl e 3. Several qua ntities of the investigated FeH lines 
dReiners & Basrill2006l) . 



surest [A] 


l°g£/ 


E, [eV] 


Br. 


J 


Q, 


magn. sen. 


9953.08 


-0.809 


0.156 


R 


10.5 


1.5 


weak 


9954.00 


-2.046 


0.199 


P 


16.5 


2.5-3.5 


strong 


9956.72 


-0.484 


0.375 


R 


22.5 


3.5 


strong 


9957.32 


-0.731 


0.194 


R 


12.5 


1.5 


weak 


9973.80 


-0.730 


0.196 


R 


12.5 


1.5 


weak 


9974.46 


-1.164 


0.108 


R 


4.5 


0.5 


weak 


9978.22 


-1.190 


0.030 


Q 


2.5 


2.5 


strong 


9979.14 


-1.411 


0.093 


R 


2.5 


0.5 


strong 


9981.46 


-1.006 


0.130 


R 


6.5 


0.5 


weak 


9982.60 


-1.322 


0.035 


Q 


3.5 


2.5 


strong 



4.3. FeH Line shapes 

To quantify the visual results of Fig. [12] and Fig. Q~3] we mea- 
sured Wa, the FWHM, and the line depth of the ten investigated 
FeH-lines (see Tab. [3). We compare 3D- and (3D)-lines with 
each other to study the effects of the velocity fields in the 3D 
models and to explore the behavior of the FeH without broad- 
ening effects from the hydrodynamical motions. A run of these 
quantities is plotted in Fig. [14] and Fig. [T31 

As we mentioned above, we have to correct the line quan- 
tities from models with changing log g for their slightly differ- 
ent r e ff. In Figs. [12] and [14] one can see how the investigated 
quantities depend on T e g. We determined spline fits p for the 
three quantities of each line. For these fitting functions p we 
took into account all five different effective temperatures. In or- 
der to correct the line quantities to a reference temperature of 

r% f - = 3275 K, we use a correction factor y = !lgl_gg } an( j 

eft ' rtjiiam.UclT-) 

multiply the quantity for the log g-model with y. This gives us 
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Fig. 12. FeH lines with constant \ogg and varying T e s. Plotted 
are pairs of 3D- and (3D)-fines shifted by a constant for bet- 
ter visibility. The models are located at T e g values of 2800 K, 
3275 K, 3380 K, 3820 K and 4000 K and a log g value of 5.0, ex- 
cept the one with T e s = 3820 K (logg = 4.9), and the one 
withr eff = 4000 K (log g = 4.5) [cgs]. 




9953.0 9953.1 9953.2 9953.3 

Wavelength [A] 

Fig. 13. FeH lines wit constant and varying log g. Plotted 
are pairs of 3D- and (3D)-lines shifted by a constant for better 
visibility. The models are located at a T e s around 3300/f and dif- 
ferent logg values of 3.0 (r eff = 3240 K), 3.5 (r eff = 3270 K), 
4.0(r eff = 3315K)and5.0(r eff = 3275 K) [cgs]. The large 
differences between 3D- and (3D)-lines stem from the hydrody- 
namical velocity fields in the 3D models. 



the value of the quantity for a logg-model which would have 
r eff = 3275K. 



4.3.1. Equivalent Width W A 

In the r e ff-series, Wa (see Fig. [T4l upper panel) decreases with 
increasing T e s . At higher T e g the number of FeH molecules de- 
crease due to dissociation and hence Wa- This can be seen in 
the 3D lines as well as in the <3D)-lines. At r eff = 2800 K, the 
influence of v.d.Waals broadening in the 3D- and (3D)-lines be- 
comes clearly visible in the line profile due to saturation of the 
FeH lines. The ten different FeH lines all behave in a similar 
manner. The only difference is the absolute value of Wa, which 
depends on the log gf- value and excitation potential Ei of each 
line. 





Fig. 14. Run of W A (top), FWHM (middle) and the line depth 
(bottom) of ten FeH-lines (see Tab. [3]) on models with different 
T e ff. Squares are the data points for the (3D)-models and dia- 
monds for the 3D-models. Wa and FWHM are on a logarithmic 
ordinate for better visibility. We connect the data points of each 
line with fitting functions (see text) to guide the eye. Plotted are 
the quantities of the 3D-models (black solid lines) and the (3D)- 
models (gray dotted lines). The models are located at T e ff values 
of 2800 K, 3275 K, 3380 K, 3820 K and 4000 K and a log g value 
of 5.0, except the one with T eff = 3820 K (log g = 4.9), and 
the one with T eff = 4000K(logg = 4.5) [cgs]. 



In the log g-series, the W A (see Fig.fTBIupper panel) increases 
with increasing logg. The change in concentration of FeH with 
lower logg, which results in smaller Wa, depends on the chang- 
ing pressure and density stratification. The difference between 
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Fig. 15. Run of W A (top), FWHM (middle) and line depth (bot- 
tom) of ten FeH-lines (see Tab. [3) on models with different 
logg. Squares are the data points for the (3D)-models and di- 
amonds for the 3D-models. We connect the data points of each 
line with fitting functions (see text) to guide the eye. Plotted are 
the quantities of the 3D-models (black solid lines) and the (3D)- 
models (gray dotted lines). All quantities are corrected to a T e g 
of 3275 K. 
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Fig. 16. 3D - (3D) corrections to FeH abundances derived from 
different FeH lines with varying quantities see Tab. [3] Upper 
panel: The models are located at a around 3300/f and differ- 
ent log g values of 3.0 (r eff = 3240 K), 3.5 (T eff = 3270K),4.0 
(r eff = 3315K)and5.0(r eff = 3275 K) [cgs]. Lower panel: 
The models are located at T eff values of 2800 K, 3275 K, 3380 K, 
3820 K and 4000 K and a logg value of 5.0, except the one with 
r eff = 3820 K (log g = 4.9), and the one with r eff = 4000 K 
(logg = 4.5) [cgs]. 

3D- and (3D)-lines at small log g-values stems from the broad- 
ening by micro-turbulent velocities and vanishes to higher log g 
values. This time the FeH lines are only mildly saturated, but the 
velocity fields in the RHD models (see Sec. 3) are strong enough 
to affect the Wa as well. Like in the r e ff-series, the ten different 
lines show no significant variations in their behavior. They only 
vary in the amount of Wa due to different log gf- values. 

Since the differences in Wa are very small, one can expect 
that the 3D correction to the FeH abundance is very small too. 
We derive abundance corrections from comparison between 3D 
and (3D) curve of growths for each set of lines on the differ- 
ent model atmospheres. The results are plotted in Fig. [16] In 
this case the correction to the different r e ff of the log g models 
is not applied. The 3D-(3D) abundance correction is between 
-0.001 dex for the coolest high \ogg model and -0.07 dex for 
the \ogg = 3.0 model. In all cases the abundance correction is 
negative which mean that the 3D lines appear stronger due to 
the enhanced opacity which becomes larger due to the micro- 
turbulent velocity. 

4.3.2. FWHM 

The dependence of the line width (measured as the width of the 
line at their half maximum (FWHM)) on T e ff is shown in the 
middle panel of Fig. [14] At low T e ff, one can see that the FWHM 
of the 3D- and (3D)-FeH lines decreases with increasing T e g. 
The v.d.Waals broadening loses influence and also the disso- 
ciation of FeH molecules leads to smaller and narrower lines. 
After r eff around 3380 K, the FWHM of the 3D lines reaches a 
flat minimum and starts to become larger again to higher r e ff. 
This rise in the line width is probably related to the rising ve- 
locity fields in the RHD models, since thermal broadening takes 
place in both, 3D- and (3D)-lines and the latter still decrease. 
The rise of the velocity in the models with a T e g of 3380 K 
and 4000 K could also be due to the slightly lower surface grav- 
ity in these models, but we think that the main influence stems 
from the higher temperatures. The (3D)-lines decrease mono- 



tonically with increasing T e f{ and reflects the behavior of Wa- 
The difference in FWHM between 3D- and (3D)-lines is very 
small at T e g = 2800 K and increases with increasing T e s to 
~ 0.8kms _1 at the highest T e g. We have seen in Sec. 3. 2 that 
this can be explained with the micro- and macro-turbulence de- 
scription. The offset between the FWHM of the ten lines is due 
to their different log gf- values i.e. large log ^/-values results in 
large FWHM. But since we are interested in broadening from 
the velocity fields, we have to take into account that lines with 
small log gf- values, i.e. weak lines, formed deeper in the atmo- 
sphere where the convective motions are stronger. These lines 
are more broadened from the hydrodynamical velocity fields and 
hence more widened. We can see in Fig.[14]that the difference in 
FWHM between the ten 3D lines becomes smaller to high T e g, 
This is also valid for the lines in the log g-series where the dif- 
ference in FWHM between the ten 3D lines becomes smaller to 
small log g- values. 

In the logg-series, the dependence of FWHM (Fig.[T5lmid- 
dle panel) is very different for 3D- and (3D)-lines. In the 3D 
case the FWHM stays almost constant with decreasing surface 
gravity between \ogg of 5.0 and 4.0 for most lines. This is due 
to the smaller amount of v.d.Waals broadening, which loses its 
influence due to lower pressure in models with smaller surface 
gravity. This is compared from the broadening due to the rising 
velocity fields. With logg smaller than 4.0, the width starts to 
increase for all lines. This increase of line width in the 3D case 
is a consequence of the hydrodynamic velocity fields which in- 
crease strongly with decreasing log g. In the (3D) case, without 
the velocity fields, the FWHM decreases with decreasing surface 
gravity and reflects the behavior of the Wa- The difference be- 
tween 3D- and (3D)-lines reach its maximal value at logg = 3.0 
[cgs] and is around 1.3kms _I . This is on the order of the ve- 
locity fields in the RHD models (see Fig. [8). One could fit ID 
spectral synthesis FeH lines to observed ones (with known T e g) 
with the micro- and macro-turbulence description (see Sec. 3.2) 
and it will be possible with the obtained velocities to determine 
a surface gravity with the help of Fig.fTTI 

We did not see any significant different behavior between the 
ten FeH lines. Neither in the r e jf-series nor in the logg-series. 

4.3.3. Line Depth 

In the bottom panel of Fig. [14] one can see the dependence of 
the line depth on T e ff. The line depth increases with decreasing 
T^, and almost all FeH lines, except the one with the lowest 
log gf values are saturated at T e s = 2800 K. The difference in 
line depth between 3D and (3D)-lines is changing in the Jeff- 
interval. At high T e ff, the line depth of the (3D)-lines is deeper 
than that of the 3D lines. At T eS = 2800 K, this difference al- 
most vanishes. This behavior is due to the saturation of the FeH 
lines at low because both, the 3D- and (3D)-lines reach their 
maximal depth. The decrease of the line depth with increasing 
r e ff is due to the dissociation of the FeH molecules at higher 
temperatures. 

The line depth of the logg-series is shown in the bottom 
panel of Fig. [15] At low log g, the line depths of the 3D and 
(3D)-lines increase almost linearly with increasing \ogg. The 
3D lines increase with a strong slope and the (3D)-lines with 
a weaker slope. The difference in line depth between 3D- and 
(3D)-models is maximal at logg = 3.0 [cgs] and vanishes al- 
most at logg = 5.0 [cgs]. It is consistent with the velocity fields 
present in the atmospheres of the RHD models broadening the 
lines and lower the line strength of the 3D lines. The (3D)-lines 
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reflect the decreasing number of FeH molecules with decreasing 
logg due to the lower pressures. 

5. Summary and Conclusion 

We investigated a set of M-star models with T e ff = 2500 K - 
4000 K and logg = 3.0-5.0 [cgs]. For these models, the 3D 
hydrodynamic radiative transfer code C0 5 B0LD was used. The 
horizontal and vertical velocity fields in the 3D models were de- 
scribed with a binning method. The convective turn-over point is 
clearly visible in the atmospheric velocity dispersion structure. 
To investigate the influence of these velocity fields on spectral 
line shapes, a description for geometrical projection and limb- 
darkening effects was applied. With the use of contribution func- 
tions, we took only these parts in the atmosphere into account 
where the lines were formed. The resulting velocity dispersions 
range from 400 - 1600 ms 1 with decreasing \ogg and with in- 
creasing T e ff from 200 - 1400 ms _1 . These values agree well 
with velocities deduced from line shapes. We expressed the hy- 
drodynamical velocity fields of the 3D models in terms of the 
classical micro- and macro-turbulent velocities. With this de- 
scription and the obtained micro- and macro-turbulent velocities, 
it is possible to reproduce 3D spectral lines on ID atmosphere 
models very accurately, hence time consuming 3D treatment of 
FeH molecular lines in the regime of cool stars is not necessary 
for line profile analysis. A comparison of our velocities with a 
set of velocities determined from observations with spectral fit- 
ting methods showed that the macro-turbulent velocities agree, 
but the micro-turbulent velocities are a factor of two or three 
smaller than the ones determined from observations. 

A line shift due to the larger up-flowing area in the convec- 
tion zone was investigated too. It is on the order of a few m/s 
up to 50 m s _I for a very low gravity model. The time dependent 
jitter in line positions is only about m/s and would be reduced 
to mm s _1 in a real star, due to high number of contributing ele- 
ments. 

In order to use FeH molecular lines for investigations of 
spectroscopic/physical properties in cool stars (e.g. Zeeman- or 
rotational broadening), we explored the behavior in a set of 
lines on \ogg and T e ff. We investigated ten FeH lines between 
9950A and 9990A on our models with the spectral synthesis 
code LINF0R3D. FeH lines react on different effective tempera- 
tures as expected due to the change in chemical composition and 
pressure. The lines showed also a weak dependence on surface 
gravity due to changing densities and pressure. The broadening 
from velocity fields in the 3D models of the log g series is very 
strong, but for the T e ^ series the broadening from velocity fields 
is almost covered by v.d.Waals broadening. The difference in 
line width for hot models is up to 0.5 km s _1 and for low gravity 
models around 1 km s _1 . That means for the ID spectral synthe- 
sis, that one has to include correct micro- and macro-turbulent 
velocities for small surface gravities or hot T e $ . Due to the fact 
that the FWHM log g dependence of FeH lines goes in the op- 
posite direction as the logg dependence of the velocity fields, 
the FeH lines become a great mean to measure surface gravities 
in cool stars. Because the velocity fields scale with logg and it 
should be easily possible to detect them. 

FeH lines with different quantum numbers do not show sig- 
nificant differences for both, \ogg- and r e ff-series. That means 
the broadening of the lines does not depend on J, Q, or the 
branch. Furthermore lines with weak magnetic sensitivity be- 
have just like lines with strong magnetic sensitivity. All lines 
are broadened in the same way by thermal and hydrodynamic al 
motions. Only the transition probability expressed in the log gf 



value influences the behavior of the lines. The line with the low- 
est ^/-values did not saturate at low T e s, but in general they are 
similar to the other FeH lines. 

It is possible to treat the FeH molecular lines with differ- 
ent quantum numbers as homogenous in the absence of mag- 
netic fields. That allows to use FeH lines to measure magnetic 
fields dReiners & Basril 120061 12007b . Hence we conclude, that 
these lines also are an appropriate mean to measure magnetic 
field strength in M-type stars. 

Appendix 




Limb darkening coefficient f 



Fig. 17. Projection factors (vertical solid line and horizontal dot- 
ted line) as a function of the limb darkening coefficient e. 

To compare the hydrodynamical velocity dispersion in 3D- 
C0 5 B0LD models with the spectroscopic quantities micro- and 
macro-turbulent velocities, we assume a simple geometrical 
model and try to resample the broadening of absorption lines. 
An intensity beam "sees" the velocity field under a certain an- 
gle and the spectral line is broadened by the projection of these 
velocities. In this sense we project the geometrical velocity com- 
ponents on a line of sight under a certain angle and integrate over 
all angles in a half sphere. We also take a linear limb-darkening 
law into account. We assume a velocity function f(v) in velocity 
space and an intensity line profile 



I(v) = V(v)-I»-L(6), 



(2) 



where V(v) is a line profile function, /'* the continuum intensity 
at the center of the disk and L(ff) a limb-darkening law. The inte- 
grated flux of the intensity line profile broadened by the velocity 
function is then given by their convolution and disk integration, 
i.e. 



F = (D f(v) * I(v) cos 6 dco, 



(3) 



with do) = sin6> d9 dip, (p e [0, 27r], 6 e [0,7r/2]. If we assume 
that V(v) does n ot vary with different positions on the disk, it can 
be factored out (Gray 2008) and 



F = V(v) * f c ■ W(v), 
with 



W( 



v) = j)f(v. 



) • L( ff) cos 6 do. 



(4) 



(5) 
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With Eq.|5]we are left with a flux-like expression for the velocity 
function. The dispersion of the velocity function is given by 



o 2 = </(v) 2 ) - </(v)> 2 

and we can write for the projected velocity dispersion E: 

22 = h I L' (<(/(v)e)2> ~ {f{v)e}2) ■ m c ° s ° dcj 



(6) 



(7) 



where /(v) 



is the velocity vector containing all ve- 



locities in a model cube, e = 



is the basis vector in 



fy(v) 
, fz(v) , 

' cos^sinfl^ 
sin (p sin ( 
cos 9 

spherical coordinates, and L(9) = 1 - e + e ■ cos(#) a linear limb- 
darkening law with the limb-darkening coefficient e. N is a nor- 
malization factor N = J 2 " d<p j^ 2 L(9) cos 6 sin 9 d9 = n( 1 - 1). 
We included the limb-darkening in the normalization because, in 
opposite to the flux, the velocity dispersion for an isotropic ve- 
locity field has to be conserved. The average is taken over the ve- 
locities and does not affect the angle dependent parts. Basically 
we have to integrate: 

((fef) - (fef = «/ 2 > - (f x f) cos 4> 2 sin 6 2 + 
(ify) ~ (fyf) sin 2 sin 9 2 + «/?) - </ z > 2 ) cos 2 + 
2«/v/v> - (fx)(fy)) cos 4> sin sin 6 2 + 
Wxfz) - (fxXfz)) cos <p sin cos 9 + 
2({f y f z ) - (fyXfz)) sin <p sin 0cos 9. (8) 

If we want to compute the average of this quantity, we need the 
mean and the squared mean of the quantities, but not the combi- 
nations of the velocity components, because these products van- 
ish in the integration over the half sphere due to their angle de- 
pendent coefficients. After the integration of Eq|7] E 2 becomes 

^ «/v 2 > - (fx) 2 + (fy) - <A.> 2 )(76 - 15) 



20(e-3) 

«/ 2 )-</ z ) 2 )(6e-30) 
20(e-3) 



(9) 



We can see, that for an isotropic velocity field f x (v) = f y (v) = 
f z (v) = f( v ) follows that cr 2 = </(v) 2 ) - </(v)> 2 and does not 
depend on e or the geometrical projection anymore. Setting e = 
0, it then follows from Eq.|9]that E 2 . = \cr\ y and E 2 = \a 2 z due 
to geometrical effects. 

The projection factors for the vertical ( 20(c-3) ) anc ' horizon- 
tal ( 20^-3) ) component are plotted as a function of the limb- 
darkening coefficient e in Fig.fTT] They vary only about 5% from 
no darkening to a full darkened disk. The reduction for the ver- 
tical velocity is about 30% and for the horizontal components 
50%. 

For the sake of completeness, we obtain in a similar way the 
mean velocities in three spatial directions {f X y, z {v)) which are 
given by 

i c 2n r"^ 2 

(fx,yjy)) projected = — I </ W (V» • L{9) COS 6 d(i>. (10) 

M Jo Jo 

The horizontal velocities vanish due to projection but there is 
still a vertical component which is reduced to geometrical and 



limb-darkening effects. 

if x( v )} projected — if y( v )} projected — 0, 

(e-4) 

</z(v)) projected = We-l) 



(ID 
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